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ABSTRACT 

In addition to its potential to probe the Inflationary cosmological paradigm, miflimeter-wave polarimetry is 
a powerful tool for studying the Milky Way galaxy's composition and magnetic field structure. Towards this 
end, presented here are Stokes /, Q, and U maps of the Galactic plane from the millimeter-wave polarimeter 
BICEP covering the Galactic longitude range 260° < £ < 340° in three atmospheric transmission windows 
centered on 100, 150, and 220 GHz. The maps sample an optical depth 1 ^ Ay < 30, and are consistent with 
previous characterizations of the Galactic millimeter-wave frequency spectrum and the large-scale magnetic 
field structure permeating the interstellar medium. Polarized emission is detected over the entire region within 
two degrees of the Galactic plane and indicates that the large-scale magnetic field is oriented parallel to the 
plane of the Galaxy. An observed trend of decreasing polarization fraction with increasing total intensity rules 
out the simplest model of a constant Galactic magnetic field throughout the Galaxy. Including Wmap data in 
the analysis, the degree-scale frequency spectrum of Galactic polarization fraction is plotted between 23 and 
220 GHz for the first time. A generally increasing trend of polarization fraction with electromagnetic frequency 
is found, which varies from 0.5%-1.5% at frequencies below 50 GHz to 2.5%-3.5% above 90 GHz. The BiCEP 
and Wmap data are fit to a two-component (synchrotron and dust) model showing that the higher frequency 
BiCEP data are necessary to tightly constrain the amplitude and spectral index of Galactic dust. Furthermore, 
the dust amplitude predicted by this two-component fit is consistent with model predictions of dust emission in 
the BiCEP bands. The polarization angles in all three bands are generally perpendicular to those measured by 
starlight polarimetry and show changes in the structure of the Galactic magnetic field on the scale of 60° . Map 
noise and systematic effects are characterized and the resulting maps and derived parameters are corrected for 
spectral mismatch leakage and time-series filtering effects. The effort to extend the capabilities of BiCEP by 
installing 220 GHz polarization band hardware into the existing BiCEP focal plane is described along with the 
results of the data analysis from the new band. 

Subject headings: Surveys — Submillimeter — cosmic microwave background polarization — cosmology: 
observations — diffuse radiation — Galaxy: structure 



1. INTRODUCTION 

Emission from the Milky Way Galaxy at millimeter wave- 
lengths is both a rich source of astrophysical information 
and a potential contaminant for cosmic microwave back- 
ground (CMB) observations. The density of gas and dust 
in the interstellar medium (ISM) varies from very low (< 
1 particle/cm^) in diffuse regions to very high (> 10^ 
particles /cm^) in molecular clouds and complexes (collec- 
tions of star-forming cloud cores at approximately the same 
distance and age). Measurements in the millimeter band have 
the potential to probe a wide range of ISM densities while 
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near-IR bands f Mart in & Whittet|1990| have enough resolu- 
tion to probe medium to high density ISM regions. In general, 
large-scale diffuse emission from the Galaxy at frequencies 
below 90 GHz is dominated by the synchrotron mechanism in 
ionized gas, and to a lesser extent free-free along with ther- 
mal emission from rotational and vibrational modes of dust. 
Emission at frequencies above 90 GHz is dominated by vibra- 
tional modes of Galactic dust ^Whittet|1992 ). The exact com- 
position and emission spectrum of the Galaxy varies across 
the sky and can be determined with multi-wavelength obser- 
vations in the infrared and millimeter measuring continuum 
emission, spectral lines, and polarization. 

Polarized radiation probes various aspects of the IS M and 
Galactic magn etic field. Optical starlight polarization ( [Davis 



& Green stein 1951 ) is due to the absorption of background 
light by dust grains aligned perpendicular to the Galactic mag- 
netic field. Polarimetric observations of stars have revealed 
that the polarization and hence the Galactic magnetic field in 
the diffuse compone nt are oriented prefe rentially parallel to 
the Galactic plane (Mathewson & Ford|[T97 ). Dust grains 
emit radiation in the infrared and millimeter bands, polarized 
along their long dimensions, making thermal dust e mission 
polarized orthogonal to the Galactic magnetic field ( [Lazar- 
ian|2007 ). While not dominant at any millimeter- wave band, 



free-free emission from electron-ion scattering can contribute 
to the measured intensity near 90 GHz; however, this emis- 
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sion is not polarized. The desire to characterize the ISM 
and the Galactic magnetic field provides the motivation for 
millimeter- wave continuum polarimetry of the Galaxy ( Hilde^ 
brand 1988; Hildebrand et al. 1999; Benoit et al. 2004; Chuss] 
et al.|2005nMatthews et al.^2 009; Culverhouse et al. 2010). 

The main goal of the Background Imager of Cosmi c Ex- 
tragalactic Polariz ation (BiCEP, Kea ting et al. (20 0la| ; [Chr 
ang et al. ( |2010| ) and its successorj^is to search for the 
umque CMB B-mode polarization pattern due to primordial 
gravitational waves, which has an amplitude determin ed by 
the energy scale of inflation ( Seljak & Zaldarriaga 1997} 
[Kami onkowski et al. 1 997| . Polarize d Galactic emission is 
an astronomical foreground ( |Bock et al.,2008) that may need 
to be confronted to make this measurement, and therefore, 
the properties o f the emission a lso moti vate the investiga- 
tion (Jo nes 2003 l|Ponthieu et al.| 2005 ; Tuc ci et al.| 2005 ; Erik- 
sen et al.i2006; Hansen et al. 2006; Amblard et al. 2007; Lar 



son etal.|2007;Eri ksen et al. 2008 ; Leach et al. 2008 ; Miville 
Deschenes eFaL|[20 08; Dodelso n et al.| |2009; Dunkley et al. 



2009). Only a few experiments have explored large-scale 
Galactic polarization properties at frequencies between 90 
GHz and 350 GHz, such as Arc heops (353 GHz , [Benoit et al. 



|2004 )), Wmap (94 GHz, Kogut et al.| (12007])), and QUaD 
(100 and 150 GHz, Culverhouse et al. (2010)). To further 
study emission from sources other than the CMB, BiCEP was 
upgraded from a two-band experiment (100 and 150 GHz) to 
a three-band polarimeter with the addition of 220 GHz capa- 
bility for the second and third seasons. This paper discusses 
unique aspects of the BiCEP Galactic observations, including 
the Gala ctic maps with the additional 220 GHz channels ( Sec 
tion|3.1| ), explains the data analysis methodology(Section |2.2 



|2.5| ), and discusses the analysis of the polarized Galactic sig 
naKSection irgp^T] ). 

2. BICEP INSTRUMENT AND MAPMAKING 
2. 1 . Brief Instrument Description 

For a complete description of the BiCEP telescope see jYoon 
let al. I (2006) and Takahashi et al. (2010). BiCEP is an on-axis 
refracting telescope with a 250 mm aperture and can scan in 
azimuth and elevation as well as rotate around the optical axis 
(boresight) of the telescope, which is less than 0.01° away 
from the center feed. BiCEP's small aperture allows the en- 
tire optical system to be cooled to cryogenic temperatures in 
a vacuum cryostat, sealed with a millimeter- wave transparent 
foam window. While located at the South Pole, the telescope 
mount is enclosed at room temperature within the observa- 
tory, protected by a fabric bellows structure. To control the 
response of the beam sidelobes and minimize ground contam- 
ination, the telescope has an inner co-moving absorptive fore- 
baffle and a fixed reflective outer ground screen, similar to that 
used in POLAR ( Keating et al. 2003b ). The readout electron- 
ics are sealed in an RF-tight cage and consist of detector AC 
biasing circuits, analog preamplifiers, lock-in amplifier cards, 
and cold JFETs in the cryostat. 

The focal plane is compri sed of 49 pairs of polarization sen- 
sitive bolometers (PSBs) ( |Jones et al.||200 2 ), two orthogonal 
detectors per feed, whose responses are summed or differ- 
enced to measure total intensity and polarization, respectively. 
For the first observing season in 2006, the focal plane had 
twenty-five feeds tuned for the 100 GHz atmospheric trans- 

list of future space, balloon, and ground CMB experiments: 
'http://cmbpol.uch icago.edu/workshops/technology2Q08/depot/meyer- | 
stephan.pdt 
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Fig. 1 . — 150 GHz FDS Model 8 dust emission prediction ( Finkbeiner et alj 
|199 9 ), shown in equatorial coordinates. BiCEP's primary CMB observing 
held is called "the southern Galactic hole", a region of low dust emission used 
for optimal B-mode detection. The two Galactic fields are used for studying 
astronomical foregrounds and Galactic physics. The Gal-weak region spans 
the Galactic plane from Galactic longitude 260° <i< 300°, while the Gal- 
bright region spans the plane from 300° < £ < 340° . 

mission window and twenty-four tuned for the 150 GHz at- 
mospheric window. For the second and third observing sea- 
sons, the focal plane consisted of twenty-five 100 GHz feeds, 
twenty-two 150 GHz feeds, and two new feeds tuned for 
the 220 GHz atmospheric transmission window (discussed in 
Appendix [A|. Three corrugated feedhom sections, cooled to 
< 4 K, couple radiation from the two high-density polyethy- 
lene lenses to each PSB while also providing a sharp low- 
frequency cutoff for the band. The high -frequency cutoff 
is defined by a set of metal mesh filters ( |Ade et aL|2006 ) 
attached to the feedhorn stack. The filters and PSBs are 
cooled to 250 mK by a "^He-^He-^He-sorption refrigerator 
system (Duband et al. 1990). Teflon filters block out-of-band 
infrared radiation to minimize optical loading. 

2.2. Scan Strategy 

Figure [T] shows the three main BiCEP observing regions 
overlaid on the Galactic dust model of IFinkbeiner et al.l 
(1999 ) (hereafter FDS) evaluated at 150 GHz. Approximately 
10,000 hours were spent observing a region predicted to have 
minimal astronomical foreground contamination in an attempt 
to detect the B-mode signature of inflationary gravitational 
waves ("CMB" region), 945 hours were dedicated to observ- 
ing the Galactic plane in a region near the center of the Galaxy 
("Gal-bright" region, 300° <£< 340°) and 1484 hours were 
spent observing the Galactic plane in a region farther from the 
Galactic center ("Gal-weak" region, 260° <i< 300°). 

BiCEP observes all regions in a similar manner. The small- 
est observing unit is a 'half-scan'; a uni-directional azimuthal 
telescope movement at constant elevation that lasts 27 sec- 
onds. To avoid potential thermal disturbances, 3.5 seconds is 
cut at the beginning and end of each half- scan. The choice 
of scan speed is bounded at low frequency by the atmosphere 
and detector stability, and bounded at high frequency by the 
bolometer time constant. Within those constraints, the tele- 
scope scan speed is chosen to be 2.8° /sec to minimize micro- 
phonics and thermal drifts. A 'scan-set' lasts approximately 
one hour and consists of 50 half-scans each in the positive 
and negative azimuthal directions at a given elevation. A cal- 
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ibration period at the beginning and end of each scan- set con- 
sists of a small elevation movement, called an 'el-nod', which 
serves as the primary relative gain calibration within a feed 
and across the focal plane. Atmospheric loading is propor- 
tional to the line-of-sight air-mass, which is well modeled as 
csc(Oei) plus an offset, where O^i is the elevation angle. El- 
nods produce a similar detector response (approximately 100 
mK peak to peak) across all detectors for an elevation change 
of one degree. To normalize the response over time, the av- 
erage response for all the detectors in each band is calculated 
and divided out. As opposed to tracking the celestial observ- 
ing center, BiCEP centers each scan-set about a fixed azimuth 
angle causing sky sources to move relative to the scan center, 
while stationary ground and scan thermal/optical contamina- 
tion remains fixed. This has the added benefit of grouping 
both the scan and ground contaminati ons into one contam- 
ination ("scan-fixed contamination"). [Chiang eFal.| ( |2010| 
takes this process a step further and removes a template of the 
scan-fixed contamination from each scan-set; however, this 
additional step was not necessary in this paper because of the 
much larger polarization signal relative to the noise. 

After each scan-set, the telescope is stepped in elevation 
by 0.25° and moved in azimuth to locate the next scan-set 
about the center of the observing region. Each set of scan- 
sets (called a "phase") consists of seven (lasting six hours) or 
ten (lasting nine hours) steps at one of four orientations about 
the boresight (0°, 180°, 135°, 315°) centered in the elevation 
range 55° to 60°. Observations of the Gal- weak field were 
carried out mostly in six hour phases in Austral winter. Ob- 
servations of the Gal-bright field consisted primarily of nine 
hour phases during Austral summer 2008, although there were 
a few other six hour phases executed at various times during 
the three years of observing. 

Timestream statistics (such as the variance, skew and kur- 
tosis of a half-scan) and el-nod calibration values are used 
to determine nominal observing conditions. These statistics 
allowed the data to be cut on various time scales such as per- 
phase, per-scan-set and per-scan bases. Gal-bright observa- 
tions use 763 out of 945 total possible hours and Gal- weak 
observations use 1463 out of 1484 total possible hours based 
on el-nod cuts, while scan statistics cut approximately 5% of 
the remaining data. 

2.3. Spectroscopic Characterization 

This section highlights the spectral characteristics of the 
BiCEP telescope, and the effects on the resulting maps. For 
a brief discussion of other BiCEP characterizations, includ- 
ing 220 GHz detector properties, see Appendix |B} for a more 
thorough discussi on of 100 and 150 GHz detectors see Taka- 
hashi et aL] ( |2010 ). The key instrumental properties of BiCEP 



are summarized in Tab le [l] as described in Appendix |B] and 
[Takahashi et aL] ( [20T0| ). 

For all PSBs used in BiCEP, at all three frequency bands, 
detected radiation from the sky produces a bolometer signal, 

d(t) = Kt^^n(t)+ g(t) J dvAe{v)F{v) j dVt P(p) 

(l+-i{Q cos(2^) + U sin(2V^)) }, (1) 

where /, 2, and U are the total intensity and two linear polar- 
ization Stokes parameters on the sky respectively. BiCEP is 
incapable of measuring the fourth Stokes parameter V, which 
accounts for circular polarization. However, for the CMB and 



TABLE 1 
Telescope Characteristics 

Characteristics of BiCEP from its three observing seasons. While not 
an exhaustive list of all possible categories and errors, the listed parameters 
show the properties of the polarimeter relevant for this paper. Section [23] 
discusses the spectroscopic characterizations in more detail, while a briet 
discussion of the other BiCEP characterizations are in Appendix |b] and Sec- 
tion l3.3.4| with further discussion in Taka hashi et al . 1 120 1 > . Ch aracteristics 
for iuu and 150 GHz are consistent with | lakahashi et aI7p|20 1 0| except for 
spectral characterizations. 



Instrument Property (Band Average) 


100 


150 


220 


Number of Feeds (2006, 2007-2008) 


25,25 


24, 22 


0,2 


Polarization Orientation Uncertainty 


<0.7° 


<0.7° 


<0.7° 


Pair-Relative Polarization Orientation 


0.1° 


0.1° 


0.1° 


Uncertainty 








Polarization Efficiency, 7 


0.92 


0.93 


0.85 


Optical efficiency (OE) 


20.8% 


19.8% 


15.8% 


Gaussian Beam Width (FWHM) 


0.93° 


0.60° 


0.42° 


Differential Pointing / Beam Size 


1.0% 


1.8% 


2.6% 


Ghost Beam Power 


0.41% 


0.50% 


1.3% 


Ghost Beam Power, Pair-Difference 


0.02% 


0.04% 


0.04% 


Spectral Band Centers, flat source (GHz)^ 


95.5 


149.8 


208.2 


Spectral Band Centers, Galaxy (GHz) 


96.3 


152.4 


212.2 


Spectral Gain Mismatch^ 


0.17% 


0.19% 


0.72% 


Relative Gain Uncertainty^ 


0.8% 


1.3% 


< 10% 


Absolute Gain Uncertainty 


2% 


2% 


15% 


NET(/iKcMBV^) 


530 


450 


1040 


NEQ per feed^ ( /xKcmb ) 


410 


340 


880 



^These values are slightly different for 100 and 150 GHz from |Takahashi et al.|(2Q10| 
due to measurement uncertainties. 

^ These values are computed before correcting for spectral gain mismatch from FTS 
measurements. 

These are the noise values used throughout this paper to calculate the white noise 
levels of the maps. 

Galaxy emission, V is expected to be negligible compared 
with the two linear Stokes parameters. A^, the effective an- 
tenna area, is assumed to be proportional to A^. To recover the 
underlying sky signal from the detector voltage timestreams 
d(t), each of the following parameters is calibrated as de- 
scribed in Takahashi et al. (20 10): the detector polarization 
angle; e, the cross-polarization response; 7 = j^, the result- 
ing polarization efficiency; PiVt), the antenna response as a 
function of angular position Vt\ F{y), the end-to-end detector 
spectral response including filters, feedhorns, lens, etc.; git), 
the detector responsivity; nit), the noise; Kt, the time-domain 
bolometer transfer function and filtering due to the electron- 
ics^ 

[Takahashi et al.| ( |2010| ) explored the leakage of total inten- 
sity to polarization for each feed, estimated using individual 
PSB pair-sum and pair-difference maps from the CMB region. 
The relative gain uncertainty is similar to spectral gain mis- 
match but can include other effects such as thermal response 
mismatch. The maps were cross-correlated, showing the rel- 
ative gain uncertainty was less than 0.8% and 1.3% for 100 
and 150 GHz feeds respectively. It was found that the relative 
gain mismatch for the 220 GHz pixels was 10% and visual 
inspection of the maps showed the two 220 GHz detectors 
gave inconsistent polarization results. The cause of this in- 
consistency was found to be mostly attributable to spectral 
gain mismatch. 

A careful campaign to characterize all of BiCEP's feeds was 
undertaken in January 2008 at the South Pole using a high res- 
olution (up to 250 MHz resolution) polarized Fourier trans- 
form spectrometer (FTS). The optical alignment and power 
falling on each pair of bolometers in a feed was calibrated be- 
fore a set of eight independent spectral measurements were 
taken. Figure |2] shows Fiy), the average spectrum for each 
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Fig. 2. — BiCEP's three electromagnetic frequency spectra, F(iy), normal- 
ized to unity (black) with the band centers (shown as vertical purple dashed 
lines). Also plotted is the average spectrum for Wmap's 94 GHz band nor- 
malized to 0.5 (gray) and the spectral radiance, S(v), for a model of atmo- 
spheric emission at the South Pole (red; Grossman (1989)), the spectrum 
(blue) of the CMB anisotropy (temperature derivative of the Planck function 
evaluated at 2.725 K ( Fixsen & Mather 2002)) and a typical Galactic emis- 
sion spectrum (green) in the BiCEP observing region, as seen through the 
atmosphere. 

band with the FTS's source spectrum divided out, assuming 
the FTS's source filled the pixel's beam. 

If not calibrated properly, mismatched spectra can cause 
spurious polarization in detector differences and gain errors 
among different fee ds. For example, the Wmap satellite 
( [Jarosik et al.|2007| notes that passband mismatch is a prob- 
lem; solved by fitting out for a spurious map component. The 
value of the mismatch is quoted as 1% on average with a 
maximum of 3.5% for the 23 GHz band (Page et al. 2007), 
in agreement with pre-flight spectral measurements. BiCEP 
does not have sufficient polarization angle coverage for each 
feed to fit out the spurious component; however, it is possible 
to use the measured spectrum to mitigate this effect. 

BiCEP's spectral response mismatch leaks intensity into po- 
larization and is caused by the combination of two effects. 
First, several feeds have mismatched spectra, most notably 
the two 220 GHz feeds. Second, each PSB is relatively cal- 
ibrated using the change in atmospheric loading with eleva- 
tion; however, atmospheric emission has a different spectrum 
than Galactic emission or the CMB. The spectral mismatch 
leakage is calculated using the measured instrumental spectral 
response, a model of the atmospheric emission at the South 



Pole ( jGrossman 1989| ), and a model of the typical Gal actic 
spectrum in the 220 GHz observing region (Section [33). 



The spectral gain mismatch (0 is : 
Ga-Gb 



Ga+G, 



(2) 



B 



where the responsivities, Ga,^, represent the two PSBs in a 
feed and are given by 



Ga,b ■ 



- atmosphere 
^Galaxy 



= / F{v)S{v)y'dv 



(3) 



(4) 



where Siy) is the spectral radiance (^ ^rm^Hz ) source 
emission spectrum of the atmosphere or typical Galactic 
source, and accounts for the throughput of the optics, 
which are assumed to be single-moded. 

Shown in Figure |3] are examples of the spectral gain mis- 
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Fig. 3. — Spectral gain mismatch for one feed at each of Bicep's three 
bands. The top plots show the measured spectrum for the two PSBs in a feed, 
"A" (Black) and "B" (Blue), along with a model of atmospheric transmission 
at the South Pole during median Austral winter conditions (red; Grossman] 
(1989)) and a typical Galactic source spectrum (green) in the 220 GHz anal- 
ysis region, as seen through the atmosphere. The bottom three plots show the 
leakage fraction per frequency, d^, if calibrated off the atmospheric emission 
but observing the Galactic source spectrum as shown. The integral of the 
leakage fraction over the bounds shown gives the total spectral leakage mea- 
sured for that feed, ^. The leakage at 220 GHz comes from the combination 
of a small difference in lower band edge and a large amount of power from 
the atmospheric water line at 185 GHz. 



match between two PSBs in a feed for each band and the 
spectral gain mismatch per frequency given by: 



dv 



= {- 



Galaxy ^B,Galaxy 



)^(^)AV(Ga + G5) (5) 



The average magnitudes of the spectral gain mismatch are 
0.17%, 0.19%, 0.72% for 100, 150, and 220 GHz, respec- 
tively, using a typical Galactic source spectrum and median 
precipitable water vapor conditions at the South Pole during 
Austral winter. Changing the observed source or atmospheric 
conditions can have an appreciable change in these numbers. 
For example, one of the 150 GHz feeds has a 0.12% leakage 
during typical Austral summer conditions but can change in 
value by 0.5% depending on the atmospheric model. In this 
paper, typical values are used to correct the maps for this ef- 
fect and simulations are carried out to probe the effects of the 
uncertainty in this calculated parameter. 

The absolute gain calibration for BiCEP maps at each fre- 
quency, co-added over all detectors, is determined by com- 
paring total intensity cross power spectra from the BiCEP 
CMB re gion to the spectra fo und by Wmap in the same re- 
gion (see [Chiang etaT] pOTOl ) and |Takahashi et al.j pOlO) for 
more details). Spectral gam mismatch between feeds could 
potentially introduce systematic effects into the analysis via 
feed calibration differences. Comparing the two 220 GHz- 
feed intensity maps, a difference of 30% is found. While this 
discrepancy was initially suspected to be due to spectral gain 
mismatch between different feeds, an investigation into the 
origin of this discrepancy did not find that this was the cause. 
While calibration per feed is important, this effect does not 
leak intensity power to polarization, and the quantities stud- 
ied in this paper are mostly relative quantities insensitive to 
this systematic. Therefore a calibration error per feed will not 
significantly affect the results of this paper. 

While the three bands are called "100 GHz", "150 GHz", 
and "220 GHz", these are not the actual band centers. The 
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band center for a given source is calculated as: 



Deg 

10 20 30 40 50 



1^0 = 



(6) 



where F{y)'v^ the measured average spectrum for each band, 
S(y) is the source emission spectrum, and accounts for the 
throughput of the receiver. For a flat spectral source, this gives 
95.5, 149.8, and 208.2 GHz for 100, 150, and 220 GHz bands 
respectively. The dominant source of uncertainty is due to the 
optical setup and whether the FTS source is beam filling or 
not. This can change these values by 0.5, 0.8, and 1.0 GHz for 
100, 150 and 220 GHz respectively. Changing the integration 
limits can also change these values by approximately 0.1 GHz 
for each band. If the band center is determined using a CMB 
source as seen through the atmosphere this gives band centers 
96.2, 150.6, and 208.8 for 100, 150 and 220 GHz. If the band 
center is determined using a typical Galactic source with dust 
and synchrotron emission, as seen through the atmosphere, 
this gives band centers 96.3, 152.4, and 212.2 for 100, 150 and 
220 GHz. Another way to calculate the average band center 
is to compute the band center for each detector separately and 
then average over all the detectors in a band. Doing this for the 
Galactic spectrum gives average band centers of 96.3, 152.3, 
212.0 GHz with uncertainties in the average of 0.1, 0.3, and 
1.0 GHz for 100, 150 and 220 GHz respectively. The standard 
deviation of the distributions over the detectors are 0.4, 1.4, 
and 2.1 GHz for 100, 150 and 220 GHz respectively, showing 
how different the spectra of a given detector in a band can 
be. Out-of-band high frequency response is less than -25 dB, 
which was characterized by using high-pass thick-grille filters 
and a chopped source. 

2.4. Time Domain Data Processing 

Following the el-nod calibration and correction for spectral 
gain mismatch, the sum and the difference of a PSB pair are 
calculated. Differencing the calibrated orthogonally linearly 
polarized detectors removes most of the unpolarized atmo- 
spheric response and unpolarized, scan-fixed, contamination. 
Atmospheric 1// noise dominates the rms of the pair-sum 
timestreams for a typical half- scan. These fluctuations can be 
less than 1 mK in good weather or reach approximately 300 
mK in bad weather, with a 1// knee as high as a few Hertz. 
The data also contain a sub-dominant scan-fixed contamina- 
tion, which does not integrate down as uncorrelated noise. To 
remove both 1// atmospheric and scan-fixed contamination, 
the data from each half-scan are high-pass filtered by remov- 
ing a second-order polynomial fit to each half-scan, at the cost 
of removing some Galactic signal as well. The pair- sum and 
pair-difference timestreams are treated separately but with the 
same filtering. 

Polynomial removal causes an obvious distortion of the sig- 
nal when the scans include the Galactic plane. To reduce this 
effect, the Galaxy is masked during the determination of the 
polynomial fit ("polynomial mask"). The preferred filtering 
scheme uses a second-order polynomial while masking out 
samples \b\ < 4°. This scheme achieves a compromise be- 
tween noise filtration and signal preservation. This process is 
not perfect and some residual filtering effects remain in the 
maps, as discussed in Section [33] A dedicated study of dif- 
ferent filtering techniques was undertaken but none improved 
the maps significantly without causing worse filtering effects 
or adding additional noise. For example, maps with only a DC 
offset removed have less filtering applied to the data but also 
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Fig. 4. — A detector timestream simulation showing three different poly- 
nomial filtering methods. The top plot shows a typical timestream resulting 
from a scan across the Galaxy (blue) with a fitted second-order polynomial 
(red) subtracted off, causing a significant distortion of the Galaxy (gray). The 
middle plot shows the same scan, except the Galaxy (green) has been ex- 
cluded from the fit. The polynomial has been interpolated across the Galaxy 
leading to minimal filtering effects. The bottom plot shows a typical scan 
that ends on the Galaxy, requiring the polynomial to be extrapolated onto the 
Galaxy. The extrapolated polynomial causes severe distortion of the maps, 
requiring these scan portions to be excluded from the analysis. 

have large scale non-physical features that makes quantitative 
analysis unreliable. 

There is an added complication in this scheme for scans that 
end within the masked region (Figure]?]). A polynomial con- 
strained by measurements on both sides of the Galactic plane 
closely approximates the low-frequency drifts within the in- 
terpolated region. However, scans not constrained on both 
sides of the plane require the polynomial to be extrapolated 
beyond the fitted region. Extrapolated polynomials tend to di- 
verge because there are no data constraining the fit. Therefore, 
the filtering scheme is modified so that scan portions that end 
within the masked region are excluded, giving rise to maps 
with a missing wedge (Figure [5]). Additionally, a brief mea- 
surement on both sides of the plane was not enough to con- 
strain the polynomials sufficiently. Therefore, scans were re- 
quired to have at least 10 samples on either side of the Galactic 
plane to be used in the mapmaking process. 

2.5. Mapmaking 

Once the data-processing steps are completed, the Stokes 
parame ters /, Q, and U are derived using standard tech- 
niques ( [Jones et al.|2007| ). Equation [T] can be simplified to 



dA^B = SA,B{l + 7A,fi (2 cos(2 + U sin(2 ) (7) 

assuming the beam functions are the same for a given pair 
of PSBs and the responsivities, gA^B^ still include the spec- 
tral gain mismatch after el-nod calibration. The Stokes pa- 
rameters /, 2, and U now represent quantities integrated over 
Vt and V, and {A, 5} refers to the two orthogonal bolometers 
within a given feed. 
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Fig. 5. — Integration time per 0.25° Healpix map pixel for 100 GHz 
(top: Celestial coordinates, second from top: Galactic coordinates), 150 GHz 
(second from bottom; Galactic Coordinates), and 220 GHz (bottom; Galactic 
Coordinates) derived from all three seasons co-added over all four boresight 
angles of Bicep observations. The wedge-shaped cuts in the map without 
integration time are fro m the omission of parts of scans with extrapolated 
polynomial fits (Section [24) . The 100 GHz celestial coordinate map is an 
example of the native observation reference frame from the South Pole. The 
red outlined area (called the "220 GHz analysis region") is an area of the sky 
where analysis at all three bands can be carried out. This region is a subset 
of the blue outlined area on the 100 and 150 GHz maps (called the "100/150 
GHz analysis region") where the final analysis can only be carried out at 100 
and 150 GHz (See Section[3). 

Before calculating the sum and difference data of the two 
PSBs in a feed, the spectral gain mismatch factors calculated 
in Equation |2] are corrected: 



dA =4x(l+0 
db = J^x(l-0, 



(8) 



where ^ is the spectral gain mismatch and the d^^ denotes 
the uncorrected timestream data. The pair- sum (d+) and pair- 
difference timestreams can then be determined: 



d,= ^^^=I+a,Q^f3,U^I 



d- = 



2 

dA - dB 



(9) 



= a-Q + p-U ^ Q cos(2 i^a) + U sin(2 V^a), 



which give rise to the Stokes parameters /, Q, and U where 
a± and /3± account for polarization angle: 



7a cos(2 tpA) ±7fi cos(2 tpB) 



7a sin(2V;A) ±7fi sin(2 V^g) 



(10) 



The A and B bolometers are assumed to be nearly perpen- 
dicular, so the pair- sum gives the total intensity to a very high 
precision. If 7 is also assumed to be equal to one, then a± and 
/3± simplify to {0, cos(2?/^a)} and {0, sin(2?/^A)}- However, 7 
is closer to 90%, and can differ by several percent between a 
given pair of bolometers within a feed, so no simplification 
is made during mapmaking using the pair-difference signal. 
After the spectral gain mismatch is corrected and after the 
sum and difference are taken, the resulting pair-sum and pair- 
differe nce ti mestreams are polynomial filtered as discussed in 
Section [231 

The /, 2, and U maps are given by: 



(11) 



which is a noise weighted linear least squares regression 
where m is the set of map pixels for a given Stokes param- 
eter, is the noise covariance matrix, and A is the pointing 
matrix. Since J(+ _) are filtered timestreams, the resulting sky 
maps Mi^i^Qiu) are also filtered. For this work, it is sufficient 
to assume that the noise is uncorrected, making N diagonal 
and simple to invert. The variance associated with samples in 
a single half-scan is assumed to be time-independent, and is 
calculated from the samples lying outside the Galactic mask 
after polynomial subtraction. For the pair-sum data, the point- 
ing matrix consists of ones and zeros, indicating whether or 
not the telescope is pointing at a particular map pixel. This 
simplifies Equation [TT] only requiring the calculation of the 
weighted mean of ^to determine the total intensity map 
mi. To determine the polarization maps mqiu from the pair- 
difference data d-, the pointing matrix consists of a combina- 
tion of the a's and /3's, and Equation [TT] can be written as: 



s 



i Cti 

il3i 



Wi 



(12) 



In Equation [T^ d- / is a single pair-difference timestream 
sample, w/ istrie inverse variance for a single sample, in- 
dex 7 corresponds to a given map pixel, k corresponds to a 
given band (100, 150, or 220 GHz), and S is the total number 
of samples from all the feeds per band in a given map pixel. 
After summing over /, the matrices per pixel per band are in- 
verted to solve for Q and IJ . For each map pixel, the seven 
quantities used to recover /, 2, and U are the weighted pair- 
sum data, weighted pair-sum hits, weighted pair-difference 
data multiplied by ol and /3, and weighted a^, and a/3. An 
eighth quantity, integration time per pixel, is also recorded in 
order to measure noise properties and observing efficiency. 
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3. RESULTS 

This section begins with a quahtative discussion of BiCEP's 
Galaxy maps, followed by a section describing uncertainties 
in the maps, a direct comparison to the maps observed by 
the Wmap satellite, and then finally a quantitative analysis of 
the properties in the maps. For quantitative calculations, two 
analysis regions were defined (Figure [5]). The "100/150 GHz 
analysis region" consists of 147 one aegree Healpix map 
pixels that have intensity values greater than zero, located at 
less than two degrees in Galactic latitude, and a polarization 
fraction magnitude less than 20%. Since there were only two 
220 GHz feeds installed in BiCEP, this limited the sky cov- 
erage for that band. Therefore, a "220 GHz analysis region", 
using 53 of the 147 pixels from the 100/150 GHz analysis re- 
gion, defines a subset of map pixels that can be analyzed at all 
three bands. 

3.1. Intensity and Polarization Maps 

Figures [6]f9] that show BiCEP Galactic maps in three differ- 
ent bands, bmned into 0.25° Healpix pixels, using second- 
order polynomial filtering while masking data \b\ < 4°. Polar- 
ization angles are defined counterclockwise from the merid- 
ian a t the map pixel, in accordance with the lAU defini- 
tion (Wler][T9^ BiCEP ab- 
solute calibration casts the maps in thermodynamic temper- 
ature, which shows CMB anisotropy with the same value at 
all frequencies. This unit is convenient for CMB analysis and 
gives a consistent reference frame for emission from sources 
other than the CMB. 

The total intensity maps in Figure [6] show the large contrast 
between the CMB temperature anisotropy (\h\ > 5°) and the 
emission near the Galactic plane. The noise in the intensity 
maps is not white; however, it is barely visible at either 100 or 
150 GHz, even off the plane, except at the edges of the map. 
The map pixels at the lower elevation have a lower signal-to- 
noise ratio due to the Healpix pixelization scheme that bins 
data into equal area pixels, while the scan strategy follows a 
Mercator projection. This effect is most readily visible in the 
220 GHz maps when comparing the upper and lower portions 
of the map. 

The total intensity maps overlaid with polarization vectors 
in Figure [7] show that the brightest Galactic emission is within 
two degrees of the plane and consists of smooth large scale 
features and compact sources along the plane. The maps show 
that the intensity signal dominates over the noise in the plane 
at all three bands and the polarization vectors are mostly per- 
pendicular to the plane. 

Figure |8] shows there is U signal in the plane correspond- 
ing to polarization vectors that are not perfectly perpendic- 
ular to the Galactic plane. There are regions of positive U 
polarization in the plane at all three bands. However, there 
is a region in the 150 GHz maps near Galactic longitude 
^ = 322° that shows a significant amount of negative U sig- 
nal that is not an artifact of the filtering or systematics. The 
noise appears mostly white over the whole observing region 
with some residual striping along the scan direction at this 
map resolution. 

Figure [9] shows there is no significant negative Q polariza- 
tion in the Galactic plane, which would have produced polar- 
ization vectors generally parallel to the Galactic plane. The 
noise is nearly identical in nature to the noise in the U maps 
because both are derived from the same pair-difference data 
with identical filtering. 
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Fig. 6. — BiCEP intensity maps from all three seasons co-added over all 
four boresight angles at 0.25° Healpix resolution, in Galactic coordinates, 
at 100, 150, and 220 GHz from top to bottom respectively. The color scale 
has been chosen to emphasize the CMB anisotropy, which is visible in all 
three bands. 

3.2. Map Quantities and Polarization Fraction Model 

T he m ap quantities studied are /, the pixel intensities (Sec- 
tion |3.5| ); q, the polarization fraction perpendicular to the 
Galactic plane; r], the exponent of a simple power law de- 
scribing the relationship between q and /(Section 3.6 ); and 0, 
the polarization angle in Galactic coordinates (Section [377] ). 
The derived variables q, u, p, and are given by: 



Q 

r 



u 

y 



1 



tan 



The main polarization quantity studied is q as opposed to p 
or w, because p suffers noise bias being a positive definite 
quantity, and there is relatively little signal in u. 

It is known that the observed polarization fraction from an 
astronomical source will be lessened by disorder in the mag- 
netic field. To explore the Galaxy's magnetic structure more, 
BiCEP maps are fit to a phenomenological power law as: 



^ = ^dmedian) ^ 



(— ) 



(14) 



where Imedian IS Calculated from the map pixels in the 220 GHz 
analysis region, r] is the slope parameter and is the 

overall polarization fraction normalization parameter. The pa- 
rameter, T], approximately represents the disorder in the mag- 
netic field which is traced by millimeter- wave polarization. 
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Fig. 7. — BiCEP intensity maps from all three seasons co-added over all four boresight angles with polarization vectors, in Galactic coordinates, at 100, 150, 
and 220 GHz from top to bottom respectively. The color scale is chosen to emphasize the emission in the Galactic plane. Polarization vectors are only displayed 
if the map pixel has a signal-to-noise ratio above 10 for intensity and above 4 for total polarization P. Polarization vectors are predominantly perpendicular to the 
Galactic plane, implying that the magnetic field in the medium sampled by BiCEP is parallel to the plane of the Galaxy. While intensity values can't physically 
be negative, the observing and analysis strategy filters the maps, causing negative values in the maps as can be seen near higher Galactic latitudes. 
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Fig. 8. — BiCEP U polarization maps from all three seasons co-added over 
all four boresight angles, in Galactic coordinates, at 100, 150, and 220 GHz, 
from top to bottom respectively. Galactic +U polarization corresponds to 
a polarization vector in the +b and +i direction. Power can be seen in all 
three bands indicating areas where the Galactic magnetic field is not ex- 
actly aligned with the Galactic plane. The significant detection of negative 
U power in 150 GHz map at Galactic Longitude 322° is physical and not an 
artifact of filtering or other systematics. 

The simplest Galactic magnetic field model to compare Bi- 
CEP data to is one where r] = 0, which has no disorder in the 
field. The measured polarization is a direct imprint of the in- 
tensity signal, related by a constant polarization factor to- 
other models describe the Galaxy's magnetic fields as be- 
ing much more random. For example, a simple toy model for 



this case is one where the Galaxy is uniformly thick and con- 
sists of a constant polarized (^o) diffuse component emitting 
with a weak intensity, /q. Scattered throughout are dense star 
forming regions with random polarization angles which, when 
integrated along the line of sight or over the beam width, will 
integrate down to very low net polarization but will contribute 
to total intensity with Ii . In this case the polarization fraction 
would be: 



Q = 



(xF 



(15) 



In this case r] = -l, implying increasing dust column density 
that contributes no additional polarized intensity. 

While both of the models here are strictly empirical, other 
studies have used similar methods. These methods involve 
fitting starlight polarization to a similar power law model; the 
exponent for polarization by absorption is related to the power 
law exponent for polarizatio n by emission in the millimeter- 
wave band by r]em = Vabs - 1 • [Fosalba et al. ( 2002 ) fit starlight 
polarization data to pys. EiB-Vy^ " ^' findi ng r]abs = 0.8, which 
implies rjem = -0.2. [Fosalba et al.|d2002] ) t hen relates thi s fit 
parameter to a magnetic field modeFfrom |Burn| ( |1966| ) us- 
ing the assumption above, namely that the polarization frac- 
tion is a product of the ratio of uniform and random compo- 
nents of the magnetic field. J ones, {19^9} fits starlight polar- 
ization data to p vs . r^"^% finding rfabs = 0.75, which implies 
Vein = -0.25. |Jones| ( |T989 ) takes the fitting process a step fur- 
ther and runs Monte Carlo simulations for observing Galactic 
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Fig. 9. — BiCEP Q polarization maps from all three seasons co-added over 
all four boresight angles, in Galactic coordinates, at 100, 150, and 220 GHz, 
from top to bottom respectively. Galactic +Q polarization (red), correspond- 
ing to a vector that is perpendicular to the Galactic plane, dominates the maps. 
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Fig. 10. — Difference between Q Gal-bright maps made at boresight rota- 
tion angles {0,315} and { 135, 180}, divided by two, in Galactic coordinates 
for 100, 150, and 200 GHz from top to bottom respectively. The left hand 
plots are the difference maps uncorrected for spectral gain mismatch while 
the right hand plots have been corrected. Some of the features in these raw 
jackknife maps are due to differences in integration time and polarization 
coverage such as near the edges of the observing area. Other features, such 
as the faint blue excess on the Galactic plane at 150 GHz map arise from scan- 
fixed or telescope systematic contamination. Most of the excess at 220 GHz 
is corrected by accounting for spectral gain mismatch while there is marginal 
change in the 100 and 150 GHz maps. There is some faint striping nearly 
orthogonal to the Galactic plane in the 150 and 220 GHz maps from residual 
1// noise leaking through the polynomial mask. 

magnetic field arrangements with varying level of randomness 
and claims a more accurate fit than using an analytic equa- 
tion. Lastly, these other works study polarization fraction 
while in this paper q is studied, which does not necessarily 
follow the same power law trend as p. For example, a Galaxy 
with completely random magnetic field directions would pre- 
dict a polarization fraction given by a power law exponent of 
T] = -0.5; however in this scenario, q would oscillate about 
zero and give an average = 0. 

3.3. Noise and Systematic Error Evaluation 

To confirm the integrity of the maps presented here, they 
were cross-checked with an independently written pipeline. 
The two pipelines produced nearly identical Healpix maps, 
reducing the possibility of coding errors or other non-physical 
errors. For example, taking the difference between the two 



pipeline's 150 GHz Q maps gave a pixel rms level five times 
smaller than the noise level, with no obvious signal features 
left. 

The data have l/f noise from atmospheric fluctuations, 
electronic readout drifts, thermal instabilities, and scan-fixed 
contamination. The polynomial subtraction removes most of 
this contamination but the polynomial masking allows some 
of this noise to leak back into the Galaxy maps. The poly- 
nomial fitting is only applied to data outside of the Galactic 
plane mask, causing not to be minimized with respect to 
the noise within the mask. However, since the noise is corre- 
lated, the polynomial fit well-approximates the noise close to 
the mask edges, but decreases in effectiveness the further the 
pixels are from the edge of the mask. Therefore, the larger the 
mask used, the worse the polynomial fit inside the mask ap- 
proximates the noise, because there is a larger gap over which 
the polynomial must be extrapolated. The smaller the mask 
used, the better the fit approximates the noise; however, this 
also increases the filtering of the signal in the plane. 

Optical imperfections, such as beam mismatch, cross- 
polarization, depolarization, and sidelobe response can cause 
systematic changes in the maps. These "telescope systemat- 
ics", studied and characterized by Takahash i et al.| ([2010 ) for 
the CMB B-mode analysis, are controlled very well partly be- 
cause of BiCEP's simple, compact design. The overall mag- 
nitude of the telescope systematics and 1// leakage in the 
Galactic maps are estimated by splitting the data into two 
halves, according to the telescope's orientation angle about 
the boresight. Angles 0° and 315° are called "boresight map 
A", while 135° and 180° "boresight map B". Since there are 
only two 220 GHz feeds, 220 GHz boresight map A is deter- 
mined from a single feed, while 220 GHz boresight map B is 
determined from the other feed. This split is the most proba- 
tive for the polarimeter because the Galactic sky coverage in 
each half comes mostly from a different set of feeds, taken at 
different times, under different weather conditions, and with 
the telescope oriented differently with respect to gravity. For 
this study the two BiCEP boresight maps per frequency band 
are studied in this paper, which can be used to gauge the gen- 
eral level of residual systematic contamination in the maps. 

3.3.1. Boresight Difference Maps 

To illustrate the efficacy of the boresight difference maps. 
Figure 10 shows the difference between the two groups for 
the Q maps for all three bands, both before and after spectral 
gain mismatch correction. This is a qualitative jackknife test, 
as it does not test quantitatively against the expected noise or 
signal leakage. Figure [T0| shows some power in the Q bore- 
sight map jackknife at all three bands, which is representa- 
tive of all the Galactic jackknife maps. The residual power 
in the 220 GHz raw maps was used to validate the spectral 
gain mismatch model and the post-corrected maps show the 
level of correction achieved. While there is still some resid- 
ual power in the 100 and 150 GHz channels, it does not affect 
the results claimed in this paper. The post-corrected 220 GHz 
maps give consistent polarization results between the two de- 
tectors; however, there is still uncertainty in the spectral gain 
mismatch parameter due mostly to changes in atmospheric 
conditions. 

3.3.2. Uncertainty due to Spectral Gain Mismatch 

Simulations were run to show the level of uncertainty in- 
duced on the quantitative parameters analyzed in this paper 
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due to uncertainty in the spectral gain mismatch calculation. 
Two sets of simulations were run; one based on the uncer- 
tainty in the atmospheric conditions and one based on the un- 
certainty in the measured spectrum. 

The uncertainty in atmospheric conditions was derived 
from computing the spectral gain mismatch using Austral 
winter and summer conditions, using a low precipitable water 
vapor (pwv) value, a median value, and a high value for dif- 
ferent observing angles. The mean and standard deviation of 
these conditions were computed, where the mean values are 
the nominal spectral gain mismatch and the standard deviation 
is used as the uncertainty. The pwv conditions used were quite 
conservative and actual observing conditions most likely had 
a much smaller range of weather conditions. The uncertainty 
from spectral measurements was calculated from the eight 
separate spectral measurements and computing the mean and 
error in the mean for each spectrum for each bolometer at each 
frequency. 

A single simulation consisted of taking the nominal spec- 
tral gain mismatch parameters for each feed and generating a 
modified list of gains using a random number generator and 
calculated uncertainty. Then usi ng th e modified list, the map- 
making steps of Section [24| and [23] are repeated, and ^(/„,,^,^„) 
and T] are found. This procedure was repeated ten times for 
each type of uncertainty, which then allowed a standard devi- 
ation from the set of simulations to be computed. 

The percent error computed from the simula- 
tions based on different atmospheric conditions is 

^^^^^ = [1.0%, 1.3%, 3.6%] and ^ = [8.9%, 15.%, > 100%] 

for 100, 150, and 220 GHz. The uncertainty from spec- 



tral measurement errors 



= [0.2%, 0.2%, 0.3%] and 



^ = [2.4%, 1.8%, 41.0%] for 100, 150, and 220 GHz. The 
uncertainty on the 220 GHz parameters are larger than the 
other two bands because the two 220 GHz feeds have the 
largest spectral gain mismatch and the smallest measured 
T] value. The spectral measurement uncertainties are five 
to ten times smaller than the uncertainties from the various 
weather conditions. The combined spectral gain mismatch 
uncertainties are relatively minor for the determination of 
^(imedianh howcvcr, thcsc uncertainties provide the largest 
single source of error on r] at all three bands. 

3.3.3. Uncertainty due to Polynomial Filtering 

Polynomial subtraction removes signal as well as noise, 
causing systematic filtering effects. To estimate the filtering 
effects, simulations are carried out to compare maps before 
and after various filtering strategies as well as to compare to 
the real data. The simulations used Wmap 94 GHz DA 1 



(Section [34]) and FDS ( ,Finkbeiner et al.| [T999) as the basis 
for map intensity. For the Wmap simulation, the map was 
first deconvolved with the beam provided by Wmap and con- 
volved with the BiCEP Gaussian beam at a given band. Then, 
the maps were downsampled from 0.125° Healpix pixels 
to 0.25° pixels. An additional simulation was carried out us- 
ing the higher resolution pixels, but this caused less than a 
5% change in the resulting simulation values. For the FDS 
simulation, tools are provided to predict the power over the 
entire sky at an arbitrary frequenc}[^ A full sky prediction 
is made at 0.016° Healpix pixel resolution for each BiCEP 
band for each measured spectral bandpass data point (approx- 

^ ^ http : // astro .berkeley . edu/ dust/index .html 



imately 200 separate frequency points per band), which can 
then be summed together to make three full sky simulated Bi- 
CEP band dust maps. The maps are convolved with the BiCEP 
beam at each band and downsampled to 0.25° Healpix res- 
olution. 

The polarization data are simulated based on these intensity 
maps by choosing parameters for Equation [14] based on Bi- 
CEP maps. For the simulations the polarization fraction mag- 
nitude was always chosen to be 2% at the median intensity 
and the power law exponent r] was varied between 0.75, 0.50, 
0.25. The polarization angle is set to 7°, the approximate aver- 
age polarization angle found in the 220 GHz analysis region, 
which sets the value of U in the simulation. Using this model, 
the polarization fraction has the possibility of taking unphys- 
ical values, therefore, the maximum polarization fraction al- 
lowed is set to 10%. A value of 25% was also implemented 
but this did not change the results appreciably. 

BiCEP pointing was used to extract /, Q, U values from the 
simulated maps to generate the simulated timestreams. The 
simulated timestreams were processed with the same filtering 
as described in Section |2.4[ leading to identical coverage as 
the real data. In addition to the nominal 4° polynomial mask, 
separate simulations were performed using 2° to 9 ° polyno- 
mial masks for the WMAP simulations. Mask values above 5° 
gave marginal improvements in recovery of studied parame- 
ters but these higher mask values also decreased the amount of 
mapped area and also created more map pixels near the edges 
of the observation area with abnormal properties. 

Polynomial masks ranging from 2° to 5° were also applied 
to the real data. Then, for each map, the same pixels were fit 
to Equation [14] and the average angle and standard error in the 
angle were computed. The conclusion from the simulations 
and comparison to the real data. Figure [TT] is that while the 
BiCEP maps have been systematically filtered, the magnitude 
of this effect is small enough not to affect the general results 
from the analysis of the maps. However, these effects do pre- 
vent additional quantitative analysis of the maps beyond what 
is performed here. The systematic filtering of the intensity is 
much better understood than the filtering in the polarization 
data. The uncertainty in the filtering correction for the inten- 
sity data is ±5%, which comes mostly from the WMAP and 
FDS inputs rather than any variation across the bands. The un- 
certainty in the polarization filtering correction is only known 
to ±25%, which comes equally from variations of simulation 
polarization parameters across each band and from the dif- 
ference between Wmap and FDS simulations. The average 
correction factors for /, ^(/,„,^,^„), r] calculated from these sim- 
ulations for the 220 GHz analysis region are 11%, 37%, 39%. 
The average correction factor for computing the average q is 
23%. These values are used to correct results for further anal- 
ysis in this paper. 

3.3.4. General Map Noise 

Background photon and detector noise dominates the white 
noise floor in the maps. The average polarization map 
pixel sensitivity can be calculated by taking the calibrated 
timestreams and computing the periodogram for each half 
scan. The average periodogram per feed for each frequency 
band is then calculated (Figure [12]). Averaging the pair- 
difference periodogram from 0.1 Hz to 1.0 Hz (correspond- 
ing to an angular size of ^ 0.5° -5°) gives Noise Equiva- 
lent Temperature (NET) per detector values of 520, 450, and 
1040 /iK,,, Vs for 100, 150, and 220 GHz respectively. After 
accounting for polarization efficiencies, these correspond to 
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Fig. 11. — Fractional change in the average map parameters vs. 
polynomial-mask size for 100 GHz (red diamonds), 150 GHz (blue trian- 
gles), 220 GHz (purple upside-down triangles), and the 150 GHz Wmap 
simulation (black circles). The parameters are intensity, Q polarization, Q 
polarization fraction at the median intensity, 77, and polarization angle. The 
units for the top four plots are relative percent difference between a given 
average map parameter, X, and the 4° polynomial-mask average map param- 
eter, X40 ; AX = {X/X40 - 1). At the bottom of each plot is the fractional 
difference, ^un filtered^ in the parameter value derived from the average val- 
ues of unfiltered simulation maps compared to that from filtered simulation 
maps with a 4° polynomial-mask. The bottom plot is the difference in the av- 
erage polarization angle compared to that derived from the map filtered with 
a 4° polynomial-mask. The simulation shows the size of filtering effects on 
the maps while the different polynomial-mask sizes show the general level of 
uncertainty in the filtering effects. The changes in derived parameters due to 
polynomial-mask size are generally smaller than the changes in the parame- 
ters between filtered and unfiltered maps. 
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Fig. 12. — Average periodogram for the pair-sum and pair-difference 
timestreams from all of the Galactic scans, after second-order polynomial re- 
moval. The pair-sum data suffers from increasing levels of 1// atmospheric 
contamination from 100 GHz to 220 GHz. Most of the 1// atmospheric noise 
is removed from the polarization data by pair differencing within a feed, re- 
sulting in nearly white noise at all three bands above 0.1 Hz. The second- 
order polynomial filtering is apparent at 0.05 Hz (the lowest frequency bin), 
where the power is lower than the white noise level. The rise at high fre- 
quency is due to the deconvolution of the detector and system time constant. 



NEQ per feed values of 410, 340, and 880 /iK,,e for 100, 
150, and 220 GHz respectively. There is very good agree ment 
between th e noise estimate from the periodograms in Taka-J 
hashi et al.| (2010 ), which used two years of third-order poly- 
nomial filtered data from the CMB region, and this calculation 
using three years of data from the Gal-bright and Gal-weak 
regions that has second-order polynomial filtering. 

The noise in a given map pixel is calculated using the NEQ 
per feed values, the integration time per 1° pixel, and assum- 
ing the integration time is split evenly between the Q and U 
maps. The average sensitivity for Q or U map pixels in the 
220 GHz analysis region is 2.8, 2.8, 16.0 /iK-rms for 100, 
150, and 220 GHz respectively. 

Another method to gauge the noise in the maps is to split the 
data in half and take the difference, which cancels off the sig- 
nal and leaves the residual noise. Maps were made from right 
going and left going scans separately, and then differenced. 
This split is used to test for detector time constant mismatch 
or general telescope thermal effects; however, it tends to be 
one of the least probative for BiCEP because the two halves 
are taken at nearly the same time, under the same weather 
conditions, and with the same set of feeds. Pixels from the 
Q and U differenced maps, with a significant amount of inte- 
gration time, were multiplied by the square root of integration 
time per pixel and the standard deviation of all the pixels was 
computed. This was done for all the pixels and the 220 GHz 
analysis region pixels. Using all the pixels, the average noise 
values are 420, 320, and 930 /iK,,e for 100, 150, and 220 
GHz respectively, which are comparable to the NEQ per feed 
noise estimate from the periodogram method indicating good 
agreement between the two methods. The noise for the 220 
GHz analysis region pixels as compared to the entire map is 
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20%, 40%, and 45% higher for 100, 150, and 220 GHz re- 
spectively, due to excess 1// noise leaking inside the mask. 

3.4. Wmap 94 GHz Band Comparison 

Wmap is a millimeter- wave satellite that has ten differenc- 
ing assemblies (DAs), all producing /, 2, and U maps over 
the whole sky using one channel at 23 GHz, one at 33 GHz, 
two at 41 GHz, two at 61 GHz, and four at 94 GHz. The com- 
bination of BiCEP data with WMAP data provides frequency 
coverage in the 220 GHz analysis region from 23 to 220 GHz. 
As an important cross check on the validity of the BiCEP ob- 
servations, a direct comparison is made between the WMAP 
94 GHz band and BiCEP 100 GHz intensity maps. 

There are two systematic differences that prevent direct 
comparison. First, the BiCEP beam is a different size than 
the Wmap beam and the maps are filtered as described in 
Section |3.2| Therefore, to make a direct comparison, each 



Wmap DA map is smoothed to Bicep's beam resolution, 
sampled into BiCEP-like timestreams, and filtered. For con- 
sistency, sampling of the WMAP maps was performed using 
both 0.0125° maps, and the nominal 0.025° maps. Since both 
sets of maps are filtered, there is no need for systematic filter- 
ing correction in this case. 

Secondly, there is a difference in spectral bandpass be- 
tween the two experiments. BiCEP is calibrated to WMAP 
by comparing the CMB fluctuations in the low astronomical 
foreground region. However, Galactic emission has a differ- 
ent spectrum than the CMB, and combined with the different 
bandpass response, this could cause a systematic difference 
between the two experiments. In a very similar manner to 
the spectral gain mismatch calculation in Section [23] the ex- 
pected miscalibration is calculated using the average BiCEP 
100 GHz spectrum, the average Wmap 94 GHz band spec- 
trum, the CMB anisotropy spectrum as seen from space and 
from the South Pole, and the typical Galactic source spectrum 
as seen from space and the South Pole. 

Different atmospheric observing conditions were used as 
in Section |2.3| The mean and standard deviation of spec- 



tral gain mismatch from the different atmospheric condi- 
tions between BiCEP 100 GHz and WMAP 94 GHz band is 
0.322% ± 0.001% (the positive sign indicates an increase in 
BiCEP power relative to Wmap). The small difference re- 
sults from the CMB and the typical Galactic source spectrum 
being very similar at this band, so even though BiCEP and 
Wmap have relatively different bandpasses, there is no spec- 
tral gain mismatch between the experiments. The extremely 
small standard deviation on this quantity results from the in- 
sensitivity of the BiCEP 100 GHz band to different atmo- 
spheric conditions, especially the emission lines outside the 
band. Instead of using the average 100 GHz spectra, this cal- 
culation was repeated for each individual detector's spectrum 
giving an average and standard deviation of anomalous gain 
factors of 0.315% ±0.055%, consistent with the first method. 

The data from the two experiments were compared in a sim- 
ilar manner to the absolute calibration routine from Chiang 
[eTal., (2010), except the comparison was done on maps in- 
stead of on angular power spectrum. The anomalous gain fac- 
tor comparing the BiCEP 100 GHz and WMAP 94 GHz band 
intensity maps is calculated as 



(MAP^MAP-l^^pBICEP^ 



220 GHz analysis region pixels weighted by the BiCEP inte- 
gration time and the WMAP and BiCEP terms represent the 
two BiCEP boresight maps and the four WMAP DAs at 94 
GHz. As a consistency check on this method, the gain was 
computed for the CMB observing region as well and com- 
pared to the absolute calibration numbers from |Chiang et al. 
( 2010| ), which found the gains using the angular power spec- 
trum, giving consistent results. The two boresight maps give 
values 5.0% and 4.7% gain increases using the 0.25° maps 
and 4.2% and 4.0% gain increases using the 0.125° maps. 
Each has a 0.5% statistical error derived from using the dif- 
ferent combinations of WMAP 94 GHz band DA maps. 

To check the uncertainty based on pixel selection or sky 
variance, a simulation was done using all pixel values \b\ < 3° 
in the BiCEP observing region (335 one degree pixels) for the 
0.125° maps. For each of 1000 different trials, 53 pixels were 
chosen at random and the anomalous gain was computed. The 
mean and standard deviation of the distribution gave 3.1 ± 
2.1% and 2.5 ± 1.9% for the two boresight maps. For com- 
parison, the same procedure was repeated for the 253 one de- 
gree pixels between -3° and -6° in Galactic latitude giving 
anomalous gain values of 1.0 ± 8.5% and -5.0 ± 6.5% 

Chiang et al. ( 2010 ) quotes the absolute gain uncertainty to 
be 2% for the BiCEP maps, which decreases the relative sig- 
nificance of the anomalous gain factor found here. The differ- 
ence in anomalous gain factor due to the bandpass differences 
was calculated to be a very small; however, this calculation 
used a simple two-component continuum Galactic spectrum. 
The real Galactic spectrum is undoubtedly more complicated, 
which could lead to a larger spectral factor difference. An- 
other likely cause of anomalous gain is from systematic un- 
certainties in t he p rocessing of raw to BiCEP-filtered WMAP 
maps (Section [T2| ), either from the beam correction or flat in- 
terpolation procedure. Irregardless of the underlying cause, 
the 4% gain difference does not represent a significant detec- 
tion of a deviation between the two experiments i n the Galaxy. 
A larg er difference was found for QUaD (Culverhouse et al. 
2010| ), despite the similarity between the instruments and 
analysis approach. However, QUaD's high-frequency cutoff 
is higher and low-frequency cutoff lower than either BiCEP or 
Wmap, making it sensitive to certain emission lines to which 
BiCEP and WMAP are insensitive. 



3.5. Intensity vs. Frequency 

By studying the spectrum of the unpolarized and polarized 
emission for a given point on the sky, it is possible to under- 
stand the composition of the ISM across the sky. The spectral 
response plots of intensity vs. frequency can show the frac- 
tion of dust and synchrotron in the Galactic plane. WMAP 
provides a vast amount of information on this topic; however, 
their highest frequency band still has an appreciable amount 
of synchrotron emission. Above 150 GHz, most Galactic 
emission comes from dust, which is where BiCEP's 150 and 
220 GHz channels add unique information. 

A common unit for distributed astrophysical millimeter- 
wave emission is differential intensity as opposed to thermo- 
dynamic temperature, where the conversion factor from ther- 
modynamic units can be calculated as: 



1 



(16) 
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where the angle brackets represent a weighted average of the and B is the Planck blackbody function. The derivative with 
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Fig. 13. — The mean intensity as a function of frequency for the 220 GHz 
analysis region for Wmap's ten DAs (gray triangles) and BiCEP's two bore- 
sight maps for each of the three bands (black squares). Bice? maps have 
been filtered by the mapmaking process, systematically lowering the inten- 
sity values. To correct fo r this, each BiCEP point and error bar has been 
increased by 11% (Section [333| . BiCEP's 100 GHz points are shifted to the 
right 5 GHz from the calculated band centers for clarity. Wmap points come 
from raw, unfiltere d maps, dow nsampled to 1 ° resolution. The noiseless FDS 
model 8 (Finkbeiner et al.|1999) predicted dust maps (green diamonds) have 
been beam smoothed, spectral bandpass filtered, and downsampled to 1 ° to 
match BiCEP's bands but not polynomial filtered. The error bars for the FDS 
points and other bands come from sky variance, not measurement errors. The 
fit to Equation [18] only using Wmap points (gray) gives a dust spectral in- 
dex Q - 5.06, larger than when BiCEP data (blue) are included, which gives 
Q = 3.70. The pink and red dashed lines are the individual positive and 
negative spectral index components from fitting to both Wmap and Bicep 
points. Bicep has two data points per band that are computed from the two 
boresight maps. The fact that two data points are nearly identical indicates 
that the systematic uncertainties are not the dominant source of uncertainty 
for the average intensity of the maps, even for the 220 GHz channels, con- 
sistent with the discussion in Section [33] Wmap points are from each DA 
separately, which can indicate the approximate level of residual systematic 
uncertainty in a given band. 



respect to temperature is: 
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Using the values of the physical constants and the BiCEP band 
centers, the conversion factors are: 



^ ^ _ ^ _ MJy/sr 
--= 0.22,0.40,0.48 ^/^, 
dT mKcMB 



(17) 



for 100, 150, and 220 GHz respectively. Figure [13] and 
the subsequent spectral analysis in this section uses BiCEP, 
Wmap, and FDS data in these differential intensity units. 

Figure[T3]shows a plot of the mean and standard error of the 
mean intensity of the analysis area pixels from both BiCEP 
and Wmap, calculated from each of the ten Wmap DAs and 
two boresight maps for BiCEP's three bands in differential 
intensity units. The points in Figure [13] are fit to 



= A, v^^ + Ad V 



(18) 



a simple two-component model. BiCEP's points and error 
bars have been increased b y 11% using the systematic filter 
correction found in Section 3.3.3] The best fit spectral index 
parameters and errors from using only the WMAP points in 
Equation [18] are [CO] = [-0.36,5.06] ± [0.02,0.33], while 
if Bicep data is included then [CO] = [-0.39,3.70] ± 
[0.07,0.11]. For reference, FDS model 8 predicts = 3.5 
in the frequency range 23 - 220 GHz. From signal simu- 
lations, the uncertainty in systematic filter correction could 
be as large as 5%; however, this factor has a relatively mi- 
nor effect on the fit. For example, a 5% uncertainty in this 
correction factor leads to an uncertainty in Q of = 1%. 
This reduced sensitivity to input uncertainty is due partially 
to the non-linear nature of the fitting model but also due to the 
Wmap points not suffering from Bicep's filtering bias. Re- 
peating the fit procedure with only one set of boresight maps, 
causes a change in Q of = 2%. Splitting the map into the 
Gal-bright and Gal-weak regions gave = [3.81,3.27] with 
cr^^ = [0. 14, 0. 15] indicating a difference in dust spectral index 
between the two regions at 2.7 sigma significance. 

Therefore, in this case, the uncertainty in this fit itself is 
larger than the uncertainty due to systematic contamination 
from weather or telescope systematics which is larger than the 
uncertainty due to systematic filtering bias uncertainty. Ana- 
lyzing the Wmap DAs separately is a way to include system- 
atics from that experiment. However, the separation between 
points in any WMAP band (See Figure 13 ) is smaller than the 
uncertainty in the average value from tnat band, making sys- 
tematics minimally important. 

The simple two-component model fits the data well; how- 
ever, additional components are not excluded. The plot shows 
the importance of measurements at higher frequency bands to 
fully understand the composition of the ISM. At 95 GHz, us- 
ing the two-component fit parameters, dust makes up 51% of 
the total emission, while at 150 GHz it makes up 87% and at 
220 GHz it makes up 97%. Finkbeiner (2003 ) gives a tem- 
plate for the full sky free-free emission based on Ha measure- 
ments. Taking their input map at 0.0625° resolution, down- 
sampled to 1 ° resolution, converted to MJy/sr using the fac- 
tors given in their table 1 at 100 GHz, the average map pixel 
value in the 220 GHz analysis region is 0.004 MJy/sr. There- 
fore, free-free is approximately 40 times smaller than either 
the dust or synchrotron component in BiCEP maps. Further 
measurements at intermediate frequencies with finer spacing 
would also provide useful information about the transition 
from the synchrotron to dust dominated regime and give more 
information about other possible emission sources. 

3.6. Polarization Fraction vs. Intensity 

The Galactic magnetic field cannot be measured directly, 
therefore measurements of the field's effect on the interven- 
ing interstellar matter are crucial. The Galactic magnetic 
field produces polarization on large and small scales across 
the electromagnetic spectrum. Millimeter-wave polarization 
measurements probe emission that spans the full extent of the 
Galactic plane making it possible to exclude or motivate mod- 
els for the Galactic magnetic field structure. 

3.6.1. Polarization Fraction Data 

One prediction of most Galactic magnetic field models is 
a trend of observed polarization fraction as a function of un- 
polarized intensity. To determine if there is a trend and what 
the nature of the trend may be. Figure [14] shows scatter plots 
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Fig. 14. — Polarized fraction, q, as a function of total intensity, /, using 147 
one degree Healpix pixels (gray) from BiCEP's three bands and two bore- 
sight maps along with Wmap 23 and 33 GHz bands. Out of the 147 pixels, 53 
were used for the BiCEP 220 GHz detectors because of the reduced observed 
region available at that band. The error bars are those due to the white noise 
in the maps and do not account for systematic errors due to the atmosphere, 
instrumental systematics, or filtering. The best fit line (blue) and constant 
model (green) are shown. The weighted average and uncertainty (black) of q 
for pixels less than the median value and greater than the median value clearly 
show the trend of decreasing polarization fraction as intensity increases. All 
BiCEP maps detect this trend with greater than 98% confidence; except Bi- 
CEP 220 GHz, which has ten times fewer feeds than BiCEP's 100 or 150 
channels. Two Wmap channels are also shown for comparison. 

of ^ VS. / for the various maps. The plotted map pixels for 
both BiCEP and Wmap are from the analysis area and the er- 
ror bars are calculated using the noise per band or DA and the 
integration time per pixel. Over-plotted on the scatter plots 
is the weighted linear least-squares fit following Equation [14] 
and the weighted mean q. Since the white noise floor value is 
used, all of the fits had reduced-x^ values much greater than 
one, indicating, to some extent, the total absolute uncertainty, 
including genuine variation about this trend across the sky. 
To better approximate the total absolute uncertainty in the re- 
sulting fit parameters, the uncertainty for each fit parameter is 

multiplied by the 
Figure p?| 



^reduced 



which are the error bars used in 



Figure 1 14] may not sufficiently support the hypothesis that 
there is a significant trend of decreasing polarization fraction 
with increasing intensity as opposed to a constant model for 
q. The difficulty is that the error in polarization fraction in- 
creases as the intensity value decreases because the uncer- 
tainty in q is proportional to , weighting the q values mea- 
sured more strongly at higher /. The fact that there is a slope 
detected for all bands from 23 GHz to 150 GHz favors a 
sloped model over a flat one (See Table [2]). 

3.6.2. Model Fit Parameters 



Figure [T5] shows a plot of the fit parameters from Equa- 
tion [T4j ana Table [2] lists these average parameters with un- 
certainties from three sources and band properties. One type 



trix, and these uncertainties have been increased by y xjeduced 
to better approximate the true parameter uncertainty and in- 
clude other noise sources such as genuine spatial variation 
across the sky. This type of error can be decreased by in- 
cluding more map pixels (possible at 100 and 150 GHz, but 
not at 220 GHz) or by increasing integration time or detector 
sensitivity. 

Second, there is uncertainty due to the dispersion of differ- 
ent maps within a band, the two boresight maps for BiCEP and 
the various DAs for WMAP. This uncertainty, especially for 
BiCEP, indicates the approximate level of systematic contam- 
ination from the atmosphere, instrument thermal properties, 
and telescope systematics other than spectral gain mismatch. 
This type of uncertainty is prevalent at all three bands and 
hinders further quantitative analysis of r]. 

Third, there is uncertainty from spectral gain mismatch that 
is quantified from simulations in Section |2.3| using different 
atmospheric conditions. This type of uncertainty compared to 
the other two is only important at 220 GHz. Overall, there 
is a strong detection of at all three bands and a sig- 

nificant detection of at 100 and 150 GHz. There is not a 
detection of r] at 220 GHz due to the combination of all three 
types of uncertainties exacerbated by the fact that BiCEP only 
has two 220 GHz feeds for only two of the three observing 
seasons (See Appendix [A| for more information). Fitting for 
T] is somewhat biased by the uncertainty in each pixel being 
proportional to /~\ as the higher intensity pixels tend to dom- 
inate the fit. For example, it is difficult to tell whether a third 
parameter is needed to fit the data properly or whether there 
are separate effects happening for lower intensity pixels that 
some models would predict. 

The top plot of Figure[T5]shows a trend of increasing ^(/,„,^,^„) 
vs. frequency with a sEarp increase between 60 and 100 
GHz. At higher frequencies where dust emission dominates 
the polarization fraction is above 2.5%, rather than below 
1.5% at lower frequency where synchrotron emission dom- 
inates. Comparing the polarization fra ction levels for th e 
Wmap maps found here to Figure 5 from |Kogut et al.|p007| ), 
which shows a histogram of polarization fraction from pix- 
els within the Galactic plane, the pixels chosen here are at 
the very low end of the distribution of polarization fractions 
within the Galactic plane. This is consistent with the choice to 
only analyze regions very close to zero Galactic latitude that 
are relatively close to the Galactic center. It's important to 
note that unpolarized emission mechanisms such as free-free 
emission contribute to the measured intensity but not the po- 
larization. This can lower the measured polarization fraction 
at bands where free-free contributes significantly, although 
that is not the case for BiCEP bands and sky coverage region. 

The bottom plot in Figure [15] shows 7/ as a function of fre- 
quency; however, no model ofthis function is proposed in this 
paper. One important point is that r] changes the relative posi- 
tions of the points in ^ vs. frequency. For example, since the 
average r] value at 220 GHz is close to zero, the value of q at 
220 GHz is relatively constant as a function of total intensity. 
Therefore, if q were evaluated at a very large value of total in- 
tensity and the plot of ^ vs. frequency remade, then it would 
appear that the polarization fraction increases relative to 150 
GHz, even though the actual polarization fraction value at 220 
GHz was nearly constant. 

The weighted mean and standard error of r] across all bands, 
assuming the BiCEP data points are corrected for filtering, is 
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Fig. 15. — ^(i,nedian) ^ as H functioii of frequency using BiCEP (red) and 
Wmap (green) data. All frequencies use the 100/150 GHz analysis region 
except for 220 GHz, which uses pixels from the smaller 220 GHz analysis 
region. The error bars are the quadrature sum of fitting errors and the uncer- 
tainty from atmospheric conditions. BiCEP's points have been corrected for 
filtering effects, increased by 37% and 39% for q and 77 respectively, based 
on signal simulations. The trend of increasing q with increasing frequency 
is apparent. BiCEP has two data points per band that are computed from the 
two boresight maps. For q^i^^^^^j the two data points are close to each other 
relative to the error bars, while for 77, the separation is close to the same size 
as the error bars. This indicates the general level of systematic uncertainties 
is not the dominant uncertainty for ^(/„,g^,a„), consistent with the discussion 
in Section [33] b ut is approaching an unacceptable level for 77. Wmap points 
are from each DA separately, which can give some idea of possible systematic 
uncertainty in a given band. 

T] = -0.47 ± 0.05. A simple, flat polarization model,/^ = 0, is 
ruled out as is the simple toy model of unpolarized embedded 
sources in a unif ormly polarized background in the Galaxy, 
from Section 3.2 7^ = -1. A more complicated Galactic mag- 
netic field model is needed to explain the measured value. A 
comparison of q and r] between the Gal-bright and Gal-weak 
regions was performed, but there was no detectable difference 
between the two regions. 

An inherent problem in analyzing the q vs. frequency plot 
is that it depends at what intensity value q is evaluated. While 
evaluating q at the median intensity gives a q value at a mod- 
erate level of intensity, other methods, such as computing the 
average of q, effectively results in finding q at the weighted 
mean intensity, as in Figure [16] The average has two advan- 
tages: a model including r] is not needed and no value of / 
needs to be chosen at which to evaluate q. For this plot, the av- 
erage from all the data is shown, neglecting different boresight 
angles and DAs. As opposed to Figure [15] this figure uses 
only map pixels from the 220 GHz region, simplifying the 
comparison of 220 GHz vs. the other bands. While the min- 
imum q occurs near 40 GHz in this plot, depending on what 
value of intensity q is evaluated, this frequency can change. 
The polarization fraction approximately quadruples from 41 
GHz to 95 GHz showing a relatively large dependence on fre- 
quency. No trend in polarization fraction above 100 GHz is 
visible, although the error bars do not exclude this possibility. 
In the end. Figure [15] demonstrates the systematics levels are 



Fig. 16. — The average and error of q from the 220 GHz analysis region 
pixels combining all the data at each BiCEP (red) and Wmap band. This 
shows a different method of computing q vs. frequency as compared to Fig- 
ure Qs] which uses a model to find the value of q, splits the data by boresight 
angleor DA, and uses a different numbero f pixels at 220 GHz. Once again a 
general trend of increasing polarization fraction vs. frequency is found with 
the minimum q occurring near 40 GHz. There is a steep increase in polar- 
ization fraction between 60 GHz and 100 GHz with the minimum fraction 
occurring near 40 GHz. 

small compared to the general trend of ^ vs. frequency as is 
evident from either Figur e [T5] or Figure [16] (even if the error 
bars are larger for Figure |16| because fewer map pixels were 
used in the average). 

3.7. Polarization Angles 

The Galactic magnetic field generally causes polarization 
angles to be perpendicular to the Galactic plane in the mil- 
limeter regime. However, BiCEP and WMAP measurements 
show the magnetic field is not exactly parallel to the plane and 
changes direction between the two observing regions (Fig- 
ure [H]). The analysis area used throughout this analysis are 
further split up into Gal-bright (27 pixels) and Gal- weak (26 
pixels) and the weighted average polarization angle for each 
BiCEP boresight maps and WMAP DA is found. The average 
polarization angle in the Gal-bright region from both BiCEP 
and Wmap is 6.2° ± 1.7° whereas the average polarization 
angle in the Gal-weak region is 20.8° ± 2.2° for both BiCEP 
and Wmap. This represents a 5 sigma detection of a differ- 
ence in polarization angle between the two regions. 

BiCEP and WMAP polarization angles can also be com- 
pared to starlight polarization data. Millimeter- wave dust po- 
larization is due to emission from particles that are pr eferen- 
tially a ligned perpendicular to the local magnetic field ( [Lazar- 
ian 2007 ). A complimentary process takes place in the opti- 
cal band due to absorption, leading to polarization angle that 
should be rotated 90° relative to the millimeter-wave polar- 
ization. |Heiles| ( [2000 | compiles a table of stellar polarization 
measurements and characteristics that can be used to compare 
polarimetry from the optical and millimeter- wave bands. Out 
of the 9286 stars in Heiles^ (1999j , only those with angle er- 
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TABLE 2 
Raw fit parameters for qvs. I 

Average fit parameters and errors for Equation[l4]using map pixels from Wmap and BiCEP in the 100/150 GHz analys is regi on (except for 220 GHz which uses 
pixels from the 220 GHz analysis region). TheBiCEP fit parameters have been corrected fo r filte ring effects (Section [333} . The systematic errors "A" for the 
BiCEP T] parameters come from atmospheric uncertainty on spectral gain mismatch (Section [2.3) . The statistical errors tor both BiCEP and Wmap come from 
the non-linear fits, while the systematic errors, "B" for BiCEP come from the two boresight maps and for Wmap, the different DAs. The median intensities are 
set by taking the median value within the analysis region (53 pixels) at a given band and rounded to two significant figures. 



Experiment Channel 


Frequency (GHz) 


^median ( MK) 


^ (^median) ^'^-^ 


Aq (%) 


V 


Aristat 




ArisysB 


Wmap K 


23 


8500 


1.5 


0.10 


-0.36 


0.05 


N/A 


N/A 


Wmap Ka 


33 


3400 


1.3 


0.10 


-0.45 


0.06 


N/A 


N/A 


Wmap Q 


41 


1900 


1.3 


0.11 


-0.48 


0.06 


0.04 


N/A 


Wmap V 


61 


880 


1.7 


0.17 


-0.56 


0.08 


0.02 


N/A 


Wmap W 


94 


800 


2.8 


0.27 


-0.60 


0.10 


0.06 


N/A 


BiCEP 100 


95.5 


800 


2.5 


0.16 


-0.64 


0.09 


0.10 


0.041 


BiCEP 150 


149.8 


1400 


3.4 


0.20 


-0.56 


0.10 


0.07 


0.063 


BiCEP 220 


208.2 


4400 


2.8 


0.40 


-0.08 


0.29 


0.23 


0.24 



rors less than 10° that are within the BiCEP observing region 
were considered. For direct comparison, the starHght polar- 
ization angles have the expected 90° difference subtracted off. 
There were 36 stars with starlight polarization measurements 
in the Gal-bright region giving a weighted average polariza- 
tion angle of -8.3° ± 3.8° while there were 24 stars in the 
Gal- weak region giving a weighted average polarization an- 
gle of 28.1° ± 2.8°. One caveat is that the starlight polariza- 
tion measurements do not exactly track BiCEP's and Wmap's 
measurements in space, as there are multiple stars in some 
map pixels and none in others. 

Since the Gal-bright and Gal-weak regions are separated 
by 60° in Galactic longitude, the change in polarization di- 
rection indicates the Galactic magnetic field has structure on 
very large scales. Further small-scale analysis is possible, but 
results are inconclusive due to the uncertainty in BiCEP's cur- 
rent polarization angle measurements. As opposed to there 
is no difference in the average polarization angle between 
pixels whose intensity is less than the median intensity and 
greater than the median intensity. No trend in polarization an- 
gle as a function of frequency is detected, showing consistent, 
independent polarization angle measurements from different 
lower frequencies where synchrotron dominates to high fre- 
quency where dust emission dominates. 

3.8. Visual Optical Depth 

In order to compare underlying astronomical objects from 
the optical to millimeter- wave band, the intensity of the emis- 
sion can be converted to a common unit used in studies of the 
interstellar medium, visual optical depth. Ay. The conversion 
can be calculated as: 



~dT 



dAy dr dB 
dr dl dT 



where the relationship between the emitted intensity and op- 
tical depth, assuming an optically thin medium, is given by: 

I = B(l-exp-^)^Br. 

A r elationship to milli mete r- wave optical de pth, r(A), derived 
by |Hildebrand| ( p^83] ) and |Dickman| ( p^78] ) is: 



Ay^ 1900 r(A) 

\250 um/ 



Assuming a dust temperature of 20 K, the conversion factor is 
calculated as: 



dAy 
~dJ 



1900 



A 



250 um 




50 100 150 

Frequency (GHz) 



200 



Fig. 17. — Polarization angle as a function of frequency for BiCEP (red 
squares) and WMAP(green triangles) for both observing regions with the av- 
erage angle for each region over-plotted (blue). Also shown are the average 
starlight polarization angles (orange) and uncertainties in the average values 
(gray shaded region) in the Gal- weak region and Gal-bright regions. The 
polarization angles do not have any trend as a function of frequency in ei- 
ther region showing consistency between the two experiments and different 
emission mechanisms. However, the angles for the Gal-bright region are sys- 
tematically lower than for the Gal-weak region. The starlight angles show a 
similar change in polarization direction between the two regions in agreement 
with BiCEP and Wmap. This rotation of polarization angle between the two 
regions indicates large scale structure in the Galactic magnetic field. There is 
good agreement between optical and millimeter-wave polarization angles in 
the Gal- weak region but poor agreement in the Gal-bright region. The BiCEP 
angles have not been corrected for filtering effects. The two BiCEP boresight 
maps per band show the general level of systematic uncertainty due to excess 
noise and telescope systematics. Wmap points are from each DA separately, 
which can be indicative of possible systematic uncertainty in a given band. 

Substituting in the values of the constants, wavelengths, and 
BiCEP band centers this gives conversion factors of: 



^^14,4.3,1.5^^. 
dT mKcMB 



(19) 



B^ 



(20 K) 



C2 (^^-1)2- 



for 100, 150, and 220 GHz respectively. One caveat to this 
calculation is that the factors are only valid for dust emis- 
sion. Section [33] showed that at 95 GHz, 51% of the emis- 
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sion comes from other sources such as synchrotron radiation. 



Using the factors of [0.51, 0.87, 0.97] from Section [33 
tensity maps from BiCEP can be converted to Ay for the dust 
component in the map as: 



dAv 



7.14,3.74,1.45 



BICEP Dust 



mK, 



(20) 



CMB 



For 100, 150 and 220 GHz. 

BiCEP detections span a thermodynamic temperature 
brightness range from [0.1, 0.3, 0.7] to [6, 8, 25] mKcMB for 
100, 150, and 220 GHz respectively. This corresponds to inte- 
grated visual optical depths Ay, using Equation|20j from [0.7, 
1.1, 1.0] to [42, 30, 36]. On average, BiCEP probes a com- 
ponent of the ISM which is more diffuse than star-forming 
regions (Ay > 10), but is more dense than the medium sam- 
pled by optical polarimetry (Ay < 2). Due to the relative 
beam sizes and dust spectrum, BiCEP measurements probe 
approximately the same density medium at local intensity 
maxima at all three bands. In principle, experiments such as 
QUaD and Archeops have the potential to probe even denser 
cloud complexes because of their smaller beam sizes; how- 
ever, Archeops had higher noise levels than BiCEP and both 
experiments have so far presented their polarization data with 
angular resolution ^ 0.5°, similar to BiCEP. 

4. DISCUSSION 

Continuum polarimetry results can be interpreted with the 
aid of a model in which the magnetic field for a given patch 
of sky is nearly constant in magnitude but has an angular 
structure which is a superposition of a uniform and a ran - 
dom component ( |Jones|1989t|Miville-Deschenes et al.|2007| ). 
The very simplest models, such as a completely uniform mag- 
netic field or a completely random one, are ruled out by Bi- 
CEP and Wmap which show a statistically significant trend 
of decreasing Q polarization fraction vs. increasing intensity 
(Section[T6|). 

The detection of dust polarization with Wmap, QUaD, 
Archeops and BiCEP constrains the degree of order in the 
Galactic magnetic field on large scales. An ordered magnetic 
field nearly parallel to the plane is detected at all millimeter- 
wave frequencies. The mean polarization angles are nearly 
constant as a function of frequency, but systematically change 
direction between BiCEP's two Galactic regions. This change 
in direction is also found in starlight polarization measure- 
ments. This change shows structure in the Galactic magnetic 
field on scales greater than 60° in Galactic longitude. 

The average degree of polarization observed in the in- 
tegrated emission from star-forming cores (Ay ^30) p < 



1.5%Stephens et al. (2010), is similar t o th e degree of po- 
larization observed with BiCEP (Section [3]6]). However, this 
coincidence does not imply that the observed BiCEP polar- 
ization arises from a superposition of unresolved star-forming 
cores, with no significant polarized component emitted by the 
diffuse medium. While there is some evidence for coher- 
ence in the magnetic field across star-forming molecular cloud 
complexes up to 100 pc in size ( Li et al.^|2009 ), as a whole, 
those complexes and cores have a nearly random distribution 
of magnetic field directions and there is no eviden ce for co- 
herence in the dense medi um on larger scales ( Glenn et al. 
1999| [Stephens et al.|2010| ). The contribution from these star- 
forming complexes to BiCEP maps would consist of polariza- 
tion with a high angular disorder, averaging to very low polar- 
ization when the beam encompasses multiple complexes. For 
example, a BiCEP 150 GHz 0.6° beam corresponds to 50 pc 



for a typical Galactic source at a distance of 5 kpc in the Gal- 
bright field, and hence poorly resolves molecular cloud com- 
plexes. Consequently, because BiCEP observes a substantial 
degree of polarization over the whole Galactic plane at ~1° 
resolution on average, BiCEP must sample a medium outside 
star-forming cores, one with an embedded magnetic field that 
retains a significant component ordered on the Galactic scale. 

On the other hand, BiCEP rules out a model where no ad- 
ditional polarization intensity comes from the higher density 
ISM. The polarization fraction drops somewhat as surface 
brightness increases, q oc /~^^^, but the decline is more shal- 
low than q oc as predicted by the toy model of Section 



3.2 



having bright unpolarized sources embedded in a polarized 
medium. Therefore, BiCEP measurements probe aligned 
grains and magnetic fields over the full range of observ ed col- 
um n densities approachi ng star-forming cores. Jones] ([1989]) 
and |Fosalba et al.| ( [2002| ) find an equivalent power law expo- 
nent from starlight data closer to r]em = -0.25, implying a more 
ordered magnetic field than what is implied from millimeter- 
wave measurements. A caveat to this comparison is those pre- 
vious works analyzed total polarization fraction p, which does 
not necessarily have the same functional relationship as q. No 
attempt to directly model q using magnetic fields is included 
in this paper. 

The detection of statistically significant polarization across 
the Galactic plane at all intensity levels shows that there are 
aligned dust grains across the entire map. The r] parameter 
indicates the nature of the Galactic magnetic field and the 
disorder within the beam and along the line of sight. The 
polarization fraction q increases as a function of frequency, 
which was somewhat unexpected; other work has indicated 
that the opposite can be true at higher latitudes due to the 
much larger polarization fraction intrinsic to the synchrotron 
emission compared to dust polarization ( Kogut et al. 2007)). 
In contrast, BiCEP has shown that dust emission can have a 
higher polarization fraction in the Galactic plane. This could 
indicate unpolarized emission mechanisms such as free-free 
play a more important role in the map regions analyzed here 
or it could be an indication that material emitting synchrotron 
radiation exists in more randomized Galactic magnetic field 
regions on average, as opposed to dust grains which may lie 
in more ordered field regions. 

The observations discussed here are consistent with a 
Galactic magnetic field whose structure and order vary in- 
versely with density. The increasing disorder in denser com- 
ponents of the ISM may be due to gravitational or dynamical 
accumulation of gas and/or feedback from star formati on. The 
high angular resolution and sensitivity of the Planck dTauber| 
et al. 20 10^) 350 GHz polarization channels should be able 
to map this structure with greater precision to discriminate 
among more detailed models for the neutral ISM. 

5. CONCLUSION 

BiCEP's Galactic observations have added new information 
and insight into Galactic physics, while also confirming pre- 
vious measurements. BiCEP samples an intermediate opti- 
cal depth of the ISM and polarization is detected everywhere 
within a two degrees of the Galactic plane with values rang- 
ing up to a few percent at 100, 150, and 220 GHz. BiCEP 
detects a significant trend of decreasing polarization fraction 
as intensity increases in the maps and increasing polariza- 
tion fraction as a function of frequency. Polarization angles 
were found to be consistent from 23 GHz to 220 GHz and 
in general agreement with polarization angles measured by 
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stars in out observing regions. Adding 220 GHz capability 
to BiCEP helped to refine the understanding of in-plane fore- 
ground emission. Improvements in BiCEP's scan strategy and 
data analysis pipeline could result in increased mapping effi- 
ciency. BiCEP data have shown astronomical foregrounds to 
be complex, and simply modeling polarization as a percent- 
age of unpolarized intensity is insufficient. While polarized 
foreground models in the Galactic plane are becoming more 
precise, less is known at higher latitudes and it is not obvi- 
ous that one can extrapolate models from the plane to these 
latitudes. Future CMB polarimeters are poised to build upon 
BiCEP's results. 
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APPENDIX 
A. BICEP 220 GHZ CHANNELS 

Prior to the second observing season, two feeds designed to observe through the 220 GHz atmospheric transmission window 
were installed into BiCEP. These feeds represent the first attempt to observe CMB and Galactic polarization in this frequency 
band using PSBs. Due to the relative sensitivity of a 220 GHz channel to thermal dust emission and its polarization over lower 
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Fig. 18. — HFSS simulation of the 220 GHz coupling feedhorn to 150 GHz bolometer. The color scale represents the electric field intensity. The phase was 
chosen to maximize the field intensity at the bolometer webs. The shape of the coupling feedhorn and distance between the last corrugation in the feedhorn to the 
position of the bolometer webs is optimized to increase the coupling efficiency and minimize the cross-polarization response. 

frequency channels, the 220 GHz feeds serve a three-fold purpose. First, to constrain the polarization and intensity of the three 
main millimeter- wave sources (synchrotron, thermal dust, and the CMB) solely from BiCEP observations, without an ad hoc 
foreground model, requires three frequency bands. Second, thermal dust emission and its polarization increase understanding 
of the physics of the Galaxy. Finally, it was unclear whether the CMB B-modes, BiCEP's primary science target, would be 
contaminated by foreground emission, and would need to be modeled and removed. The 220 GHz feeds can carry more weight, 
per feed, than the other two bands for measurement of dust polarization due to the steep spectral index of millimeter- wave dust 
emission. More 220 GHz feeds were planned to be added in 2008 but logistical constraints at the South Pole in 2008 ultimately 
prevented this, limiting the experiment to using two feeds for the two final observing seasons. 

Emission from the Galaxy has the potential to contaminate the faint B-mode signal in all observing fields. Since dust emission 
follows a thermal spectrum, by monitoring polarized interstellar dust at higher frequency than the primary B-mode bands (100 
and 150 GHz), it is possible to monitor dust contamination with higher signal-to-noise per feed. While BiCEP's CMB region, 
located at high Galactic latitude, certainly exhibits minimal dust column-density and therefore low-intensity millimet er- wave 
emission, this doesn't imply it will exhibit low-polarized emission due to the details of magnetic grain alignment (Lazarian] 
[2007 ). Little is known about the polarization properties of thermal dust emission in the high-frequency bands above 100 GHz 
used for CMB observations. BiCEP's 220 GHz feeds are a unique link between the low-frequency CMB polarimetry of Wmap 
and the high-frequency survey of the Archeops experiment; from centimeter wavelengths to the sub-millimeter. The combination 
of Wmap, QUaD, Bicep, and Archeops polarization observations of the Galactic plane cover more than a decade of frequency. 

BiCEP PSBs were only fabricated for 100 and 150 GHz operation. The 220 GHz pixels therefore were modified from 150 
GHz PSBs using in the first season of observing. The 220 GHz feed installation consisted of replacing two 150 GHz feedhorn 
stacks with a set of 220 GHz feedhorn stacks. Each feedhorn stack consists of three separate feedhorn pieces, all of which used 
the same outer- forms as the 150 GHz feedhorns. Two of the feedhorn inner-forms were made with a similar profile to the 150 
GHz feedhorns. The most difficult part of the design process was to optimize the feedhorn that couples the radiation to the 
bolometer. In this case, the 220 GHz coupling feedhorn and corrugation profile had to be matched to the existing 150 GHz 
bolometer and housing. Using HFSS[^ the feedhorn-to-bolometer housing coupling, co-polar beam shape, total throughput, and 
cross-polarization beam response were optimized. Specifically, a straight 220 GHz feedhorn profile inserted into the 150 GHz 
housing resulted in the best overall simulated performance (Figure [l^. The lower frequency cutoff of the 220 GHz band is 
defined by the feedh orn waveguide c utoff while the higher cutoff is aefined by a set of metal mesh filters, borrowed from the 
ACBAR experiment ( |Runyan et al.|2003| ). 



B. TELESCOPE CHARACTERIZATION 
Summarized here ar e the characterizations of parameters in Table [T] not mentioned previously. For more information, see 



Takahashi et al. ( 2010| ). The detector polarization angle, ip, is measurea with a 0.7° uncertain ty across all three bands, while the 
relative angle uncertainty in the pair is measured with a 0.1° accuracy ( ^Takahashi et al.|20T0) ). Nominally this is the mechanical 



orientation of the bolometer web in the focal plane, however in practice, this is a parameter characterized after deployment for 
more precise analysis. The ability to measure this parameter with such low statistical and systematic error is one of BiCEP's 
design advantages. 

The average cross-polarization response, e, measures the relative magnitude of the residual signal in the orthogonal or cross 
polarization direction given a purely polarized input in the co-polar direction. The HFSS simulations predicted the 220 GHz 
feeds should have an e less than 0.07, similar to the measured value. Depolarization is larger for the 220 GHz channel due to the 



Ansoft's HFSS: http://www.ansoft.com/products/hf/hfss/ 
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coupling effects between the 220 GHz coupling feedhorn and the bolometer and housing optimized for 150 GHz observations. 
The resulting polarization efficiency, 7, is the loss in polarization signal compared to the unpolarized signal used for calibration. 
7 directly affects the observed polarization fraction, however, it has been measured to 2% uncertainty, which is much smaller 
than other errors in the measurement. 

The optical efficiencies (OEs) were derived by taking load curves while observing blackbody radiation at different temperatures 
and give the total end-to-end sensitivity on the sky. While this parameter is not used in the mapmaking pipeline, it is implicitly 
included in the telescope noise estimates. The OEs for the last observing season were 20.8%, 19.8%, and 15.8% for 100, 150, 
and 220 GHz respectively. HESS simulations for the 220 GHz feeds predicted a coupling efficiency over 98% (similar to the 100 
and 150 GHz feeds) with 2% reflected power. The OE of the feeds that were converted from 150 GHz to 220 GHz declined by 
an average of 15%, while the other feeds did not change between seasons. The simulated transmission and reflectivity did not 
include the lenses, which were optimized for 125 GHz operation (midway between 100 and 150 GHz), not 220 GHz. 

The beam functions, P(ft), were measured using a circularly polarized broadband noise source for 100 and 150 GHz feeds 
and Jupiter for 220 GHz feeds. The beam response functions are fit with an elliptical Gaussian giving an average beam Eull 
Width Half Max (EWHM) of 0.93°, 0.60°, and 0.42° for 100, 150, and 220 GHz respectively. The differential beam size and 
ellipticities are very small for BiCEP and do not affect the results of this paper. Both polarized and unpolarized sidelobe response 
are characterized by observations of a broadband source and are determined to be negligible for this analysis. 

The absolute telescope pointing is calculated from observatio ns of m ultiple stars taken by an optical camera mounted to the 
telescope. The absolute pointing is the same as used in Chiang ^ al.| ( [2010) and is independent of a given feeds transmission 
band. Relative radio pointing for each feed is calibrated using individual maps of CMB temperature anisotropy. The differential 
pointing mismatch between PSB pairs has similar values across all three bands, with a mean of 0.0041°, 0.0046°, and 0.0047° 
(1.0%, 1.8%, 2.6% of the beam size) for 100, 150, and 220 GHz respectively. The 220 GHz intensity maps show some small 
evidence for differential pointing between the two PSBs within a feed and between the two feeds. 

Reflections in the optical system created a "ghost image" opposite to the primary image with respect to the boresight for each 
feed. The average ghost image power, as determined by observations of the Moon, was 0.41%, 0.50%, and 1.3% for 100, 150, 
and 220 GHz respectively relative to the primary image. The magnitudes of the ghost beam images relative to the primary image 
are nearly identical for a given pair of PSBs in a feed, giving power differences of 0.02%, 0.04%, 0.04% for 100, 150, and 220 
GHz bands respectively. Since the Moon's high intensity lowers the responsivity of the detectors, the ghost image systematics 
given here are upper limits. These differences are below the noise level in the Galactic polarization maps for all three bands. 
While this effect has been measured well using a bright source such as the Moon, there is no evidence for ghost effects in the 
Galactic maps. 

A miscalibration of the time-domain impulse response, Kt, would act like a beam size mismatch leaking intensity to polariza- 
tion. However, the time-domain response for each of the three bands is measured to a very high precision and deconvolved from 
the raw 50-Hz-sampled timestream. There is no evidence for time constant mismatch from either the Galaxy data or CMB data. 



